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■ Abstract 

^| , Influence of the nuclear binding in the deuteron deep inelastic structure function 

F2(x,Q 2 ) is investigated. The description is based on the Bethe-Salpeter formalism in 

■ the ladder approximation and the operator product expansion within an effective meson- 
nucleon theory. It is shown that the binding in a fully relativistic treatment can be simply 
parametrised and its origin can be understood in terms of conventional nuclear degrees of 

■ freedom. A numerical estimate of the EMC effect in the deuteron is given. 

t-H ■ 

^ ■ 1 Introduction 

Ever since nuclear effects in deep inelastic lepton scattering on hadron targets have clearly 
manifested itself over the last decade, both particle and nuclear theorists obtained a profound 
insight into physics of the EMC effect. However a majority of models of the EMC effect provided 
by a fully relativistic hadronic theory resorts to the Relativistic impulse approximation (RIA) 
^ \ with the off-mass-shell kinematics, for example see Refs. |], |2[. 

It seems to be evident that the free nucleon structure function (SF) F^(x,Q 2 ) is somehow 
affected by the nuclear medium. Here Q 2 and x are the four momentum transfer squared and 
the Bjorken scaling variable x = Q 2 /{2mv) ) with v being the energy transfer in the laboratory 
frame. Conventional nuclear physics models, for a review see Ref. 0, assume that properties 
of nuclear constituents, nucleons and mesons, are insensitive to the nuclear environment. Yet 
the hadronic current, which is poorly understood object so far, may depend on dynamics of 
nucleons inside a nucleus and their interactions. The common assumption in this case is to treat 
the hadronic current as the sum of currents of single (quasi)free nucleons, whereas the relation 
between the momentum and energy of the virtual nucleon is specified in a relevant treatment of 
the relativistic bound state. Such a proviso substantiates derivations of the convolution formula 
for the nuclear SF with a consideration of relativistic kinematics. 

The approach we take in this paper follows a rigorous method introduced in Refs. 0, |5| to 
investigate deep inelastic scattering on light nuclei. The method relays on the leading order 
operator product expansion (OPE) within the effective meson-nucleon theory. It allows one 
to derive the convolution formula for the nuclear SF which accounts for the role of the Fermi 
motion and binding effects in a new and unambiguous form. Such a formulation yields the 
nth moment of the nuclear SF as a product of the moments of the isolated-hadron SF's and 
the moments of the effective hadron momentum distributions. The latter are given by matrix 
elements of twist-two operators sandwiched between the ground state vector of the target. 

■"■On leave from Laboratory of Theoretical and Nuclear Physics, Far Eastern State University, Vladivostok, 
690000, Russia 
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Using the Hamiltonian formalism relevant operators in the OPE are found in the lowest 
order approximation on the meson-nucleon coupling constant and resulting operators involve 
the meson degrees of freedom. Finally the matrix elements of these operators are calculated 
over the ground state of the target in the non-relativistic (NR) limit. The coefficient functions 
together with meson cloud contributions are substituted by the moments of the isolated-hadron 
SF. 

In this paper we apply the method proposed in Refs. ||, |[ in order to describe the deuteron 
SF F®(x, Q 2 ) where the system is viewed as the relativistic bound state in the Bethe-Salpeter 
(BS) formalism. To this end expressions for the twist-two operators in the OPE are found 
by making use of the relativistic Hamiltonian. Explicit expressions for the moments of the 
deuteron SF are found through actual calculation of the matrix elements in terms of the BS 
vertex function. At the first step we assume that the off-shell effect on quark and anti-quark 
momentum distributions in hadrons may be neglected^]. In this way the Mellin transform 
reconstructs the deuteron SF as a sum of convolution terms describing the Fermi motion and 
the binding in the deuteron due to meson exchanges. Since meson exchange currents (MEC) are 
associated with exchange mechanisms responsible for the binding in the deuteron, momentum 
distribution for the constituent mesons reveal itself as well. 

Advantage of this relativistic approach is that it differentiates the pure binding effects from 
the off-shell effects as both have an unlike origin. The magnitude of the former is given in 
terms of well-understood quantities such as the kinetic energy of a nucleon and the binding 
energy of the deuteron. Besides dealing with the deuteron we can obtain all required formulae 
analytically within well-defined approximations and subject the problem to intense scrutiny. 

The paper is divided into five sections. In Sec. II the basic formalism is reviewed. The twist- 
two operators in the OPE of the Compton amplitude is specified. In Sec. Ill the moments of 
the fractional momentum distributions of nucleons and mesons in the deuteron are calculated. 
The NR limit of the energy-momentum sum rule is performed. Computation of F^(x,Q 2 ) as 
a function of x is carried out in Sec. IV. And Sec. V contains a summary of principle results 
and conclusions. 

2 Structure functions of the deuteron in the BS ap- 
proach 

2.1 The relativistic vertex function of the deuteron 

Consider a system of two nucleons interacting via a scalar, pseudo-scalar and vector meson 
exchanges. The Lagrangian density for meson-nucleon couplings is 

d = -gvN^NV - g s NNS - ig^N^N^ + ... (1) 

with N being the nucleon and S, (pi and meson fields. 

Assuming that the NN interaction to proceed through the meson exchange only, we start 
with the one-boson-exchange (OBE) model in application to the BS equation (BSE) in the 
bound state region of the deuteron. Therefore the BS amplitude, xip'i P) = (®\TN Pl N P2 \P), 

2 Though it may well be incorporated in this line of considerations. 
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and equivalently the BS vertex function F(p; P) = S 1 (pi,P2)x(p] P)i satisfy the BSE in the 
momentum space as 

rfa P) = i J 0^VobeAp - pW; P), (2) 

where S(pi,p2) is the free two-nucleon Green's function and Vobep is the OBE interaction 
kernel. The BSE (|j) is written in terms of the relative momenta p, p' of the nucleons and the 
center-of-mass (cm.) momentum P and p\$ = P/2 ± p. The interaction kernel Vobep is the 
sum of one-particle exchange amplitudes of certain bosons with given masses and couplings 
generated by Ci in Eq. ([[]). 

Since the BSE is homogeneous, it cannot determine a multiplicative constant of x(p] P). In 
order to normalise the BS amplitude, condition based on the electro-magnetic (EM) current 
conservation at zero momentum transfer can be used. In the IA the normalisation condition 
reads 

-jjj d A p X (p;P)[ Wo S- 1 (pi,P2)]x(p;P) = Po, (3) 

where P = VP 2 + M 2 and \ stands for the BS amplitude conjugate of x- An appropriate 
normalisation constant in Eq. @ is denoted as N. In the Euclidean space the amplitude x an d 
its conjugate x are related as follows: x(p<i, p;P) = — [x(~ PA: P! P)]*1^1o'\ where p 4 lies on 
the imaginary axis. 

The Lagrangian (|l|) and the BSE (0) define an effective meson-nucleon theory for the de- 
scription of the bound state of the deuteron which is meaningful in a low-order approximation. 
Such a theory is not fundamental one in a sense that interactions in Eq. ([I]) are viewed as 
effective ones which are modified by vertex form factors. If the BSE is applied in the ladder 
approximation as done in Eq. (fj), all further calculation of matrix elements of EM currents 
expressed in terms of the nucleon and meson fields over the deuteron state have to be made up 
to the second order in the strong coupling. This procedure is justified in describing the nuclear 
structure effects in electron scattering. 

2.2 Operator product expansion 

The structure functions can be calculated via the OPE of the amplitude for forward 
Compton scattering of virtual photons off a hadronic target |], [7j]. The expansion separates 
characteristic hadronic and nuclear scales and formulates this particular analysis of the nuclear 
SF as a calculation of matrix elements of certain operators within low-energy physics. 

The spin-averaged Compton amplitude is expanded in terms of local operators O and coef- 
ficient functions E in the form 

T^(P,q) oc - g, v £ {^-] g w ...g Mn £ El a (Q 2 , g)(P\0^\P) 

n=even \V / a=N,B 




ff w E El a (Q 2 ,9)(P\0^\P), (4) 

a=N,B 
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where O^ 1 " ^" are series of operators having twist two at the level of naive dimension counting 
and a tags the fundamental fields of the theory, i.e. nucleons (N) and bosons (B). Each operator 
is appropriately symmetrised and trace-subtracted. In this procedure the deuteron SF F 2 is 
expressed in terms of invariant diagonal matrix elements fx^ D 

(P\0^-^\P) = {P m . ..P^}fi a J D , (5) 

and the coefficient functions E 2 a by 

M n - 1 (F 2 D )= £ E®(Cf,g)rit D , (6) 

a=N,B 

where M n (F) defines the nth moment of F, M n (F) = Jq dxx n ^ 1 F(x, Q 2 ). In Eq. (|6]) we neglect 
the effects of nucleon and target masses corrections which can be taken into account in a closed 
analytic form. The explicit line of our calculation proceeds as follows: (1) Define a basis of the 
twist-two operators in OPE (Q). (2) Parametrise the coefficient functions in terms of moments 
of SF's F 2 of nucleons and bosons, F^ a oc M n _ 1 (i 7 ^ 1 ). According to our assumption E™ are 
identified with the corresponding moments of a measured isolated-hadron SF, which in case of 
a = N include effects of meson clouds surrounding a nucleon. (3) Obtain the moments of the 
deuteron SF by calculating the matrix elements \i^[ D in Eq. (El). (4) Use the Mellin transform 
to invert the moments and exhibit SF's. 

Let us determine explicitly the operators in the expansion (^) for the meson-nucleon theory. 
For a nucleon field and a pseudo-scalar (and scalar) field interacting via coupling given by (HD 
there are only two series of symmetric bilinear operators 

-J S{: iV 7 m d ■ ■ ■ d N :}, (7) 

/ ^ \ n— 1 ^ n 

0^"=[-) S{:<pid ...d (8) 

where S implies complete symmetrization in the Lorentz indices fi\. . . fi n and removes all traces 
with respect to each pair of indices. 

If a theory involves fundamental massive vector bosons, matters become more complicated 
by the fact that the vector field itself has twist zero and some subsidiary condition is needed 
to introduce the set of relevant operators. Two sets of the operators can be formed. 

Set (a): 

-j S{:N^D ■ ■ ■ D N:}, (9) 

= U S{:F x ^d ...d F/ n :}, (10) 

l 
2 

and set (b): 



Cy'"^ = ( i ] ^{: V^ 1 a ... a :}, (li; 



0W~^ = (_j S{:N<f*d ...9 iV:}, (12) 



= ( _ j S{:V x d ... d V x :}. (13) 
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Set (a) is introduced on the ground of gauge-invariance ||. The first two operators labelled 
by iV and F are invariant under the transformation V^(rr) — > V^(x) + gvd fi A(x) and N(x) — > 
exp (igyA(x))N(x), although the mass term in the free Lagrangian for the vector field — \\XyV 2 
breaks this gauge symmetry; and the subsidiary condition is used that the symmetrised 
form of the energy-momentum tensor for the massive vector bosons is restored by the operators 
when n = 2: 

qi» = F \i* Fx « + ^ V VV V - ^g v {NYNV u + N^NV") - g^C. (14) 

The operators Oy 1 "' fln ([T3J) guarantee that. On the other hand set (b) is chosen on the basis of 
that the Lagrangian of the vector field is taken in the form of a covariant sum of four Lagrangians 
each separately corresponding to the components of and imposing the invariant subsidiary 
condition d ■ A = 0. Thus the operators depending explicitly on the interaction and composed 
out of the tensor are excluded from the set. In this case the energy- momentum tensor is 
restored in the form 

W = d^A x d v A x - g^C (15) 
and the operators (O) and (13) preserve the momentum sum rule for the F 2 (x), i.e. M 1 (F 2 ) = 1. 



3 The moments of the deuteron structure function 

3.1 Nucleon contribution 

The set (a) was considered in Ref. in a calculation of the nuclear SF within the Walecka 
model. The result is that the exchange by vector meson is incompatible with the observed 
magnitude of the EMC effect. The nucleon part of the nuclear SF, because of the final-state 
interaction with the vector meson field of a nucleus, is found to be twice as large as the 
experimentally observed one. Moreover the meson part peaks at small values of x, which is 
either not observed. 

We consider the twist-two operators that belong to the set (b) . In this set operators for the 
vector and scalar meson fields are taken into account on equal basis. Taking the Bjorken limit, 
defined by Q 2 \v — > oo with x fixed, and choosing a coordinate system so that q = (z/, 0j_, <#), 
where q\ = —\/Q 2 + v 2 , leads to the following. The mathematical part is to compute the nth 
<9_ derivatives, q ■ d — z/<9_, in Eqs. (0), @ and (|i~2| ) and (|T3| ) by making use of the equation of 
motion for the coupled meson and nucleon fields. Consequently the resulting operators involve 
the meson degrees of freedom. The matrix elements of the output operators over the deuteron 
state can be calculated by the Mandelstam method UlO |. 



A representation of the moments M n (F^) in terms of diagrams is shown on the Fig. 1, 
where the inner solid lines denote the nucleon propagators and the wavy lines denote the 
meson propagators. The closed triangle and square indicates the corresponding moments of 
the isolated-hadron SF. The moments of effective nucleon distributions, ^ L \ are described by 
Fig. 1(a) and (b) and can be written down as 

u^ = u^(RIA) + u^(int), (16) 
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(a) (b) (c) 



Figure 1: The moments of the deuteron SF up to the second order in meson-nucleon 
coupling: (a) relativistic impulse approximation with a struck nucleon projected to its 
mass-shell, (b) and (c) meson exchange currents. Closed triangles and square denote 
the moments of the free nucleon and meson SF's, respectively. The vertical dashed 
lines separate operators and the BS vertex function of the deuteron in corresponding 
matrix elements. 



where we call the first term in Eq. flT6|) as the RIA with the on-mass-shell kinematics^ 

/if D (RIA) = — / d 4 px{p;P)l^S-\p 2 )x(p;P)i^) n -\ (17) 

where S~ l (p2) = m — p 2 and the +- component of a vector defined as the sum of the time and 
longitudinal components originates from p ■ q m up + ; p is on-mass-shell momentum, p 2 = m? . 
The form of Eq. (|17) is close to that of Ref. ||11|| . Although the nucleon of the upper-half circle 
at Figs. 1(a) and 1(b) is confined kinematically to the mass-shell, its propagation is described 
in a covariant way as given in Eq. (|17|). 

The contribution /if' (kit) in Eq. ( |TED is due to the strong interaction between two nucleons 
in the deuteron and is described by Fig. 1(b). It is obtained through the calculation procedure 
of the operators (§) in the Hamiltonian approach: 

/if D (int) = ^ (18) 

1 r A^A^' r> / P++P+~ fc + \n-l 
1 r dpdp „ ; 01(1) V „ (V\ v(v J. P \P± ~ ^ 2 ) 



1 r A^nA^nl fi' n ~ l ( P'++P++ k + \n-l 

i- ^^^P)VoBEP(^)[7V]'"x(p-;f) ^i^^^ ; ■ 



P + J (2ti)' '" ™ ™> ' ^-^-1^ 

where k = p' — p and p, p' are the on-mass-shell nucleon momenta. 

The first two moments of [i^ D have a physical interpretation. For n = 1 we find that 
/if D (int) = and as a result the baryon number associated with the nucleons fields is preserved. 
The normalisation condition fl3|) implies that /if D = /ii(RIA) Ar / D = 1. For n = 2 Eq. ( |T8| ) 
contributes to the energy-momentum sum rule. The fraction of the total momentum of the 
deuteron carried by nucleons is 

^> = ^(RIA)+^ /J, (int), (19) 



3 Below the normalisation constant N is included in the definition of the BS amplitude. 
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where at the deuteron rest frame, P+ = M, we have 

^' d (ria) = ~ Jd l P p + x(p; Ph+s-'btixlp; n (20) 

, r(int) .__L / to P) ^ w + wvp^ + P) . (21) 

In order to estimate contribution given by Eq. (|21|) we can manipulate with the BSE (0) 
and find 

/V ( mt ) = Jfi J d PX(P;P) g (Pa)x(P;-P)- ( 22 ) 

The nucleon contribution to the deuteron momentum can be written as (2) = 1 — A, where 
A is the missing fraction of the deuteron momentum carried by exchange mesons. Notice that 
Eq. (p2] ) suggests that the magnitude of the binding in the deuteron does not explicitly depend 
on the interaction kernel. 



3.2 Meson exchange currents 

The meson part of the nuclear SF is concentrated at small values of x. We expect that the 
bulk of the relevant meson effects in deep inelastic scattering on the deuteron to be associated 
with pion exchange. Therefore, we consider the role of pions only in this section. Although 
other kinds of meson exchanges which appear in the OBE kernel of Eq. (|2]) are taken into 
account in this approach, it is impossible to estimate their contribution to deuteron SF. There 
is no data on the free SF's for these mesons. We can also expect that contribution due to the 
scalar and vector exchange in the deuteron cancels, if both kinds are considered in the manner 



given by the operators in Eqs. (§) and (|13|) . 

In the impulse approximation, the contribution of the mesons according to Fig. 1(c) is 
summed up to (here we write down expression for 7r-mesons only): 

k+ , k< 



u n/D 



J ^kn^k;P)^±(^r-\ (23) 



where Vl\ = k 2 — fi 2 and function n lx {k) can be interpreted as the excess pion 4- momentum 
distribution associated with the one-pion exchange in the deuteron 

n ^=^j^^W x(p+t ' p) ' (24) 

where the one-pion exchange kernel V n appearing in Eq. (U) is used in solving the BSE @. 

The first two moments of yJlJ D are interpreted as the mean number of the excess pions in 
deep inelastic scattering, jj^ D = (N n ), and the missing fraction of the deuteron momentum 
carried by the excess pions, = (y}^, respectively. Since the twist-two operators from 

the set (b) for n = 2 are related to the energy-momentum tensor, the deuteron SF obeys the 
momentum sum rule [0, [J: 

1 

F°(x)dx = ^(P|e 00 + e 33 |p) = 1, (25) 
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where the components of the energy-momentum tensor correspond to those in Eq. flTol). 
Eq. ( p5|) ensures that the MEC carry away the missing part of the deuteron momentum. 

3.3 Non-relativistic limit 

In order to find the NR of the theory, though it may seem to be dubious, we need to take into 
account only the positive-energy states in the BS amplitude. The negative-energy components 
do not survive in the limit |p>[/m. <C 1. 

We write down the BS amplitudes in the rest frame of the deuteron as decomposition in 
terms of the product of Dirac spinor of nucleons: 

xip; P (0) ) = «(p)«(-p)[^(p)w(po, IpI) + ^JUpMpo, IpI)], (26) 

where ii(p) is the positive-energy Dirac spinor. We rewrite Eq. (|26|) equivalently as 

X (p; P (0 )) = «(p)u(-p)*(po, p), (27) 

x(p;P(o)) = -w(p)s(-p)**(-po,p), 

where the radial wave functions u(po, |p|) and w(po, |p|) can be expressed in terms of the BS 
vertex functions T L with L = 0, 2 as follows (E p = Vp 2 + m' 2 ) 

i I i\ i r (o, |p|)^o(po, |p|) . . h 1 r 2 (o, |p|)6(p , IpI) , oq , 

«(Pd,|p|) = ^ [(f-^ p ) 2 -^] ' ^'IPl) = ^ [(M.^.pg] ( 28 ) 

with a smooth function £l(po, |p|) reflecting the dependence on the relative energy p of T^. 
Let us consider Eq. fll7|) for the positive-energy components only. We find that 

tf' D QUA) = ijj d 4 P Mp)\ 2 (l + ^){E P - pl){^T-\ (29) 

where p+ = (E p , p) and p^ = (Ma/2 — p , — p). Our next step is to integrate over the relative 
energy picking up the pole in the 4-quarter: 



dpo £ L (Po, IpDMpo, IpI) £l(p)£l>(p) 



Po =E p -f 



(30) 



^ ,2vr (f-E p -p + ze)(f - £ p + Po + z e )2 (M _ Ep - p ) 2 
Keeping in mind Eq. ( |3~0|) we reduce Eq. ( f2~9"D to the form 0] 

^ /D (IA) = / dpn(p)(l + ^)(^ i ) n " 1 (31) 

with the nucleon density n(p) = |\1/| 2 depending on radial wave functions of the L = 0, 2 angular 
momentum states (£ ~ 1) 

^ r i( o,| P |) 



iV VM(2E P - M) 
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which in the NR limit are usually identified with the conventional wave functions in 3 Si- and 
3 .Di-states. Eq. ([H]) is the impulse approximation with a hit nucleon on its mass-shell. 

The other important quantity is (z), because its magnitude is well-known from NR studies. 
First of all we integrate over p in Eq. ( p2|) in the same manner as done above. We get 

ftf D (int) = ±J dpn(p)(M - 2E P ) = A( M - 2(E p )). (33) 

Using Eq. (|3T| ) for n = 2 we obtain 

^ /D (RIA) = A«£ p > + i^). (34) 

It is clear that the proton and neutron in the deuteron carry the fraction of the total 
momentum of the deuteron given by |J 

IA int 

where (T) = (p 2 )/m is the kinetic energy, (V) = — (T) + Ed is the potential energy of two 
nucleons in the deuteron and Ed = M — 2m. 



4 Results and discussions 

Owing to Eq. (§) the SF F 2 D (x,Q 2 ) is written as a sum of convolutions of the SF of the 
unbound nucleon and pion and the momentum distribution functions of the nucleons and pions 
inside the deuteron (we assume that the SF F^ ,7T do not carry any dependence on the relative 
4-momentum). The contribution of the interaction corrections to the RIA given by Eq. (|18|) 
can be approximated in a very plain form yielding 

F 2 D (x, Q 2 ) = F 2 N/D (x, Q 2 ) - { ^J-F 2 N (x, Q 2 ) + 6F^(x, Q% (36) 

where (V) = m/i/y D (int) and 

F 2 N/D (x,Q 2 )= J f N / D (z)F 2 N (x/z,Q 2 )dz, (37) 

z>x 

6F^ D (x,Q 2 )= J r/ D (y)FZ(x/y,Q 2 )dy. (38) 



y>x 



The distribution functions f N ^ D (z) and f n ^ D (y) describe the momentum distributions of the 
on-mass-shell nucleons and off-mass- shell pions carrying the fractional longitudinal momenta 
of the deuteron y and z, respectively. We have at the rest frame: 

f N,D (*) = ^J d4 P P (o))! 2 ( X + |r) 5 { z ~ §) > ( 39 ) 
F /D (y) = Jd'kn^P^sL-^Y (40) 
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where p a = (pE p , p) is the momentum of the struck nucleon in the partial state a = 2S+1 L P J=1 ; 
4> a denotes a partial amplitude. 

Both functions f N ^ D (z) and f n ^ D (y) specify the number densities of the nucleons and excess 
pions and satisfy the normalisation integrals 

f N ' D {z)dz = 2, (N n ) = f r ID {y)dy. (41) 



The SF F2 probes the momentum distribution of the constituents and thus conservation of the 
energy-momentum requires that 

f' D {z)zdz + / r /D {y)ydy + [ f B/D (y)yd y = 1, (42) 



where f B ^ D (y) the momentum distributions of some other mesons in Vobep- 

For the sake of a simple analytical analysis of the nucleon density in Eq. fl39|) and leaving 
no doubts of accuracy of theoretical predictions, we employ the BS vertex function which is 
the solution of the homogeneous BSE for two spin-1 particles interacting through a covariant, 
separable potential Hl2| . The relevance of such an approximation to the expressions with the 



original BS vertex function of the OBE potential could be justified. Indeed the final formulae 
for /i^ D (int), (z) and (Nb) (v. i.) do not contain the OBE potential. On the other hand 
discarding of the negative-energy components in the vertex function is also justified by the 



conclusion in Ref. 12 



The magnitude of the binding in the deuteron can be calculated in a straightforward way 



by making use of Eq. (22). Computation of (z) in the energy-momentum sum rule gives: 
the on-mass-shell nucleons carry (z)ria = 1-0153 which breaks down the sum rule (pE2"p. The 
interaction term ([18]) takes away a fraction of the longitudinal momentum, /z 2 (int) = —0.0188, 
and reduces (z)ria yielding the final value, (z) = 0.9965. 

Unfortunately the total contribution of the mesons cannot be calculated explicitly without 
running into formidable computation. Yet we find by means of Eq. (|42| ) that the fraction of 
the momentum carried by the exchange mesons is A = 0.0035. The first moment ^ D is the 
mean number of the exchange mesons associated with the reaction. The exact expression being 
similar to Eq. (||) can be approximated by the formula 

(Nb) ^^rj d 4 px(p;P)—S- 1 (p 1 ,p 2 ) X (p;P), (43) 

where p_ = p — p z . It is expected that (N B ) is totally accounted for by the pions, since 
the part due to heavier meson exchanges is suppressed in average by their larger masses in 
the propagator VL 2 k and their effective couplings. If we pretend not to take interest in a detail 
structure of the f 1T ^ D (y) that emerges from Eq. (f40"D, we can tentatively estimate the momentum 
distribution of the pions for not sufficiently large values of z by exploiting results of Ref. [13| . 



The quantitative analysis shows that results are not sensitive to the precise functional shape of 
f n / D (y) but controlled mainly by (N^) and (y). Thus f n ' D (y) could be parametrised imposing 
the balance between (N^) and (y)^ in the simplest manner. 
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Figure 2: The ratio of the deuteron and 
isoscalar nucleon structure functions. Solid 
curves: curve 1 denotes the relativistic im- 
pulse approximation (RIA) with on-mass- 
shell kinematics; curve 2 denotes the RIA 
with taking into account the nuclear binding; 
curve 3 denotes the sum contribution of the 
RIA, the binding and meson exchange cur- 
rents. The dashed curve corresponds to the 
solid curve 3 with a different parametrisation 
of F^(x,Q 2 ). Experimental data are taken 
from Ref. j\L$. The error bars are combined 
statistical and systematic errors. 
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4.1 Comparison with data 

Now we are in a position to perform numerical analysis. Fig. 2 shows the EMC effect 
in the deuteron at Q 2 = 25 GeV 2 given by the ratio R(x,Q 2 ) = F 2 D (x, Q 2 )/F^(x, Q 2 ) in 
the present approach. All computations of the nucleon momentum distribution f N l D (z) and 
/i2(int) are carried out using the BS vertex function in the framework with multi-rank separable 



interaction 12. The free nucleon = (F 2 P + F 2 n )/2 and pion F£(x } Q « 25 GeV 2 ) SF's 



obtained from the combined proton and deuteron data of the BCDMS collaboration |15| and 



data on massive lepton pair production of the CERN NA10 collaboration [16|], respectively, are 
taken from Refs. || and [|13[ . 

At small x, x < 0.2, the MEC increase the deuteron SF relative to the nucleon one. At 
exactly x = the ratio has maximum, R(0) = 1 + (N n ) F£ (Q) / (0) , where (N n ) » 0.015, and 
falls off roughly linearly as x grows. The value of F£(0)/F^(0) is known from the parametri- 
sation of the SF's. As can be seen, the nucleon contribution F N ^ D (x,Q 2 ) to R(x,Q 2 ) is not 
important in this range of x. 

At intermediate x, 0.2 < x < 0.6, the MEC are negligible as long as f 1T ^ D {y) vanishes at 
y — > 1, but the correction due to strong interaction between nucleons, v. s. Eq. (|36|) , becomes 
responsible for the dip. This depletion is totally controlled by the value of /^(int). In a more 
realistic case we still may have contribution from deep inelastic scattering on pions in this 
region, because the argument x/y of F% in the convolution formula, see Eq. ([38]), will run over 
broader range in f n ^ D (y)- 

The rise of the ratio R(x) at large x, x > 0.6, is associated with the relativistic Fermi 
motion. The tendency is that the harder "tail" of F®{x) in the vicinity of the boundary for 
the single-nucleon kinematics, the harder the momentum distribution f N l D (z) in Eq. (ft7|). 

Unfortunately pertinent data on the ratio R(x, Q 2 ) are not still available. However in 
Ref. |L7j it is suggested that the A-dependence of the F^(x) should be determined by the local 
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properties of nuclear matter, i. e. the average nuclear density p(A). The naive extrapolation to 
the deuteron case leads to R(x, Q 2 ) <=a 1 + [F^jF® — I) / {p{A) / p{D) — 1), although this (may) 
overestimate the effect due to the isoscalarity of the deuteron. The model-dependent value 
of the ratio F^ / estimated from averaging over measured data on various nuclei F^/F^ 
in SLAC-E139 experiment is presented on Fig. 2. Within this hypothesis the deuteron has 
a significant EMC effect which contradicts to our theoretical predictions together with the 
parametrisation of the neutron SF ||: a deviation of R(x, Q 2 ) from unity at intermediate x is 
extremely large, about 4 % versus the theoretically predicted 2 %, and centers of data points 
continue to decrease up to x = 0.7 which is the feature of data for heavy nuclei. 

Our computation analysis shows that R(x, Q 2 ) does not exhibit substantial dependence on 
Q 2 in the broad range and on a specific parametrisation of the neutron SF (both effects are of 
the second order relative to the binding). The ratio is more sensitive to the behaviour of the 
free nucleon SF as x — ► 1, since this region is enhanced by the nucleon momentum distribution. 



5 Concluding remarks 

In this paper the deuteron SF F2(x,Q 2 ) is examined in the theoretical approach based on 
a picture in which nucleons and mesons are the relevant degrees of freedom. The deuteron is 
treated as two spin-| particle bound state through a consideration of the Bethe-Salpeter equa- 
tion in the ladder approximation. In the model the deuteron SF given by Eq. (|36|) is obtained 
as a sum of three terms. The first term known as the relativistic impulse approximation is the 
convolution of the the free nucleon SF with the fractional momentum distribution function of 
the nucleons projected to their mass-shell. The second term in Eq. ( p6|) accounts for the binding 
correction to the relativistic impulse approximation. Similarly to the first term the last piece is 
the convolution of the free pion SF with the excess pion fractional momentum distribution in 
the deuteron. The momentum distributions of the deuteron constituents are derived from the 
matrix elements of certain twist-two operators in the operator product expansion method and 
expressed in terms of the BS vertex function. 

The principal result of the paper is that the binding in the deuteron in the approach in- 
corporating the relativistic dynamics in terms of the conventional degrees of freedom could be 
revealed in an unambiguous way. The relativistic impulse approximation is obtained in the 
form that avoids the explicit dependence of the momentum distribution on the relative energy 
of nucleons and thus we are able to present accurate estimates in a manner similar to the 
non-relativistic treatment of Ref . || . 

The magnitude of the nuclear binding and Fermi motion correction together with a neutron 
SF F£ shape x-dependence of the EMC ratio for the deuteron. High precision measurements of 
the nucleon SF at large x combined with a reliable theoretical treatment of the deuteron may 
help to resolve the question of the extraction of the neutron SF from experimental data on a 
deuterium target and a heavy nuclei target 



The open questions the model suggests are that an additional work to be done in order to 
establish connection between the off-shell and on-shell SF F^ in this framework; next to what 
extent the negative-energy states in the BS vertex function impact on the momentum distri- 
butions of the nucleons; and finally the exact functional shape of the momentum distribution 
of the excess pions associated with the one pion exchange in the deuteron is not of the least 
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interest. 

To sum up, we can conclude that the point of view developed in this paper might help to 
discriminate the contribution of the pure binding effects in the nuclear SF and regard corrections 
due to the relativistic Fermi motion in a transparent way. 
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